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the  vibrational  relaxation  process  of  a  molecule  plays  an 
important  role  in  the  processes  of  high  speed  gasdynamics, 
chemodynamics  and  ultrasound,  especially  in  the  molecular  laser 
til.  A  lot  of  relaxation  experiment  data  have  been  accumulated 
£21.  However,  only  some  macroscopic  physical  quantities  are  able 
to  be^  measured  in  experiments;  or  instance,  in  a  shock  tube  or 
laser  florescence  experiment,  the  relaxation  process  of  a 
corresponding  vibrational  mode  can  be  monitored  by  observing  the 

variation  of  a  spectral  radiation  with  £ime  [3,  41.  This 

|  eu'^Xl  ■'■^iit 

variation  follows  the  form  of,  I  -  I  =  I^e  O  ,  except  in  the 


beginning  short  moment.  The  decay  Constants,  x  >  contain 
information  of  various ' relaxation  rates  of  the  system.  However, 
this  does  not  simply  represent  the  rate  of  one  particular  process, 
except  for  a  pure  biatomic  molecule.  This  paper  presents  the 
topic  of  how  to  determine  the  relaxation  rate  from  decay 
constants. 
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Simplified  processes  are  adopted  by  most  of  the  past  works. 
The  models  used  in  analyses  are  too  imperfect,  and  application 
range  is  restricted;  moreover,  Schafer  formula  [51  which  is  used 
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to  process  V-T  data  has  no  relation  with  the  one  used  to  process 
V-V  data. 


A  universal  data  process  method  is  presented  in  this  paper, 
and  a  formula  which  is  suitable  for  various  concrete  conditions 
(such  as  single  A  and  double  k  )  is  thus  obtained  .  The  methods 
for  both  V-V  and  V-T  rate  processes  are  also  unified.  This 
explains  that  the  Schafer  formula  is  an  approximation  of  this 
universal  formula  under  special  conditions.  Inaccurate  data 
presented  in  Ref.  [2^.43  are  also  pointed  out  by  giving  an 

example  of  CO  system. 
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1.  Method  and  Universal  Formula 


First,  for  an  experiment  system  a  relaxation  equation  set 
describing  various  energy  transferring  processes  is  listed  based 
on  a  certain  relaxation  mechanism.  Then  these  equations  are 
expanded  with  respect  to  a  relative  bias, 
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A  relationship  between  the  exponent  constant,  A,  and  the  rate,  K,  _ _ 

r 

is  thus  obtained  by  solving  the  characteristic  equation,  f(A,K  ) 

*  □ 

of  the  linearized  equation  set.  The  function  of  K  =  K  _  C  A,  □ 

i  i 

fC  )  can  then  be  solved  reversely.  The  desired  K  can  be  _ 

J  *i  i 

calculated  by  a  given  K  and  the  experiment  data,  - 

i  '  - 

f  Cod«s 
Avtiil  and/or 
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Based  on  a  CO  system  as  an  example,  an  equation  set  is 

2 

given.  In  CO  -N  -H  0  mixture,  its  relaxation  mechanism 
2  2  2 

is  *r,  *„<— »CC*(>0.  co/iy—  *,  +  n),  COjC^,*— *2i*»),  CO/*,)  <—*7 

Its  relaxation  equation  set  is 


where  the  relaxation  terms  are  defined  in  Ref.  Cl,  63.  The 
associated  solution  of  the  linearized  equation  set  of  X  is 


C3) 


where  C  is  a  constant  associated  with  the  initial  conditions, 
ij 

The  characteristic  equation  which  satisfies  the  exponent  constant 
is  a  fourth  order  polynomial  of  X  . 


Since  a  Fermi  resonance  exists,  the  «*-*2»*  rate  is  high. 
Moreover,  because  precise  and  complete  data  for  this  process  is 
lacking,  it  is  practical  to  choose  an  approximation  in  which  *t 
and  %  are  in  balance.  In  other  words,  the  rate  is  much 

higher  than  X  and  other  rates.  Therefore,  the  characteristic 
equation  becomes 
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where 


In  process  A, 
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where  f  Is  gram  molecule  fraction,  #  is  vibrational  characteristic 
(V) 

temperature,  N  is  the  particle  number  in  unit  volume. 

Equation  (6)  is  the  summation  of  all  components  of  the  system,  and 
K  is  the  relaxation  rate  of  this  process. 
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Under  a  resonance  condition,  it  is  not  difficult  to  prove  that 
the  spectral  radiation  intensity,  I,  of  which  the  jumping  quantum 


number  is  1,  is  proportional  to  the  vibrational  energy  E  of  that 

particular  mode  C73.  In  the  experiment,  the  time  history  of  the 

spectral  radiation  density  is  measured,  and  its  decay  constant  is 

then  determined.  Obviously,  this  decay  constant  should  perfectly 

correspond  to  the  eigen  values  of  X  for  that  particular  mode. 

i 

If  some  of  the  eigen  values  are  large,  their  effect  on  Equation 
C3)  will  be  reduced  to  zero  in  a  very  short  period.  Thus,  for 
those  systems  containing  multiple  eigen  values,  only  one  of  two 
decay  constants  can  be  measured  from  experiments.  For  a  single 
ljthe  desired  rate  can  be  calculated  from  Equation  C4)  by  giving 
this  f  and  other  known  rates,  and  then  the  unknown  KtA;  N1  can  be 
determined  from  Equation  C6)  by  a  giving  S  and  K.tA;  Cl  which  is 

obtained  in  a  pure  CO  experiment.  If  twoU,  can  be  measured  in 

2  '  ' 

an  experiment,  then  both  f(X  )  =  0  and  fCX  )  =0  have  to  be 

1  2 

solved  simultaneously  to  obtain  two  desired  rates. 


In  a  pure  CO  and  a  CO  -H  0  systems, 

2  2  2 

on  lyuM****  appear  in  the  relaxation 

process.  Take  the  characteristic  equation  when  *  and  n  are 

1  2 

in  equilibrium,  which  is  equivalent  to  take  „  I  —  9  in 

Equation  (4).  For  a  single  , 


If  y|-*3»}is  chosen  as  a  V-V  exchange  mechanism  in  a  CO 

2 

molecule,  + *1  in  Equation  (2)  has  to  be  substituted  by 

Compare  the  approximation  of  the  characteristic  equation  which 

assumes  if  and  #'  in  equilibrium,  and  with  Equation  C4)  and 
1  2 

(7),  there  exists  a  simple  relationship  as 
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2.  V-T  Rate  Formula  in  a  Limiting  Case 


The  Schafer  Formula  of  V-T  relaxation  time,  *•*,  at  the  minimum 
vibrational  energy  level  of  a  system  which  is  converted  from  the 
total. relaxation  time,  measured  in  some  conventional 

experiments  Csuch  as  ultrasound  and  a  shock  tube)  151 
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where  Cv  is  the  total  vibrational  constant  volume  specific  heat; 
CvC  ^)  is  the  part  contributed  by  mode.  This  formula  is  stil 
being  extensively  used. 


Actually,  the  V-T  rate  formula  can  be  derived  from  the 
limiting  case  when  the  V-V  rate  is  much  greater  than  the  V-T 

rate.  For  the  pure  CO  case,  the  final  form  of  Equation  (7)  is 

2  l  l 

similar  to  Equation  (9)  in  the  case  of  - >—  After 

fw*i  r» 

substituting  the  expression  form  of  the  vibrational  specific  heat 


m 

n 

Cr>>| 

ft 


1| 


m 


into  Equation  C9),  it  clearly  indicates  that  the  Schafer  formula 

is  an  approximation  form  in  this  limiting  case.  From  Table  1,  it 

shows  that  this  approximation  is  good  for  CO  when  the 

2 

temperature  is  low. 

Similar  conclusions  can  be  drawn  for  other  systems  and 
relaxation  mechanisms. 


Table  1  Comparison  of  the  Schafer  formula  and  the 
limiting  case  of  Equation  C7) 
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0.3  **4 

9.33M 

0.6472 

0.60*12 

0.7007 

0.6691 

0.9440 

0.9194 

Key:  (l)  Factor*  of  Equation  (7);  C2)  Factors  of  Schafer; 
(3)  Relative  error  <%). 


3.  Results 


Table  2  and  3  show  the  results  of  Case  K _  . «  •  K  and 

J  •  TB 

fC  are  obtained  from  the  approximation  of  Ref.  C2,  31  and  Ref. 
RWG 

C43,  respectively,  or  extracted  directly  from  literature.  K.  is 
the  result  of  this  paper  (data  of  decay  constants  are  taken  from 
Ref.  [3,  41). 


It  shows  in  these  tables  that  the  differences  can  be  as  large 

as  40%,  except  in  the  CO  -H  0  case.  Since  errors  caused  by 

2  2 

the  temperature  variation  are  different,  additional  bias  can  be 
introduced  into  those  temperature  dependent  relations. 
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4.  Discussion 


The  base  of  the  method  presented  in  this  paper  is  to  linearize 

a  nonlinear  relaxation  equation  set.  Both  experimental 

measurement  and  theoretical  analysis  indicate  that  X  decays 

i 

exponentially  with  respect  to  time.  This  implies  that  the  second 
order  higher  terms  approach  to  zero  much  faster  than  the  first 
order  term.  This  also  means  that  effects  of  the  high  order  terms 
on  the  signal  are  noticeable  only  at  the  beginning  short  moment; 
however,  what  we  are  interested  is  the  decay  situation  after  a 
long  period  of  time,  and  this  should  mainly  depend  upon  the  linear 
term.  The  estimation  which  is  drawn  from  the  conditions  of  an 
laser  fluorescence  experiment  indicates  that  the  linear  term  still 
plays  the  major  role  even  in  the  initial  highly  non-equilibrium 
stage.  This  linearization  approximation  is  especially  good  for 
those  experiments  which  have  small  bias  from  equilibrium  (X  << 


Both  TB  12,  33  and  RWG  [43  are  processed  for  single  X  and 
double  fc  ,  respectively.  Let's  analyze  their  approximation 
characteristics.  It  is  not  difficult  to  prove  that  the  formulae 
used  by  them  can  be  derived  from  the  original  equation  presented 

in  this  paper  and  takes  ! ►  oo  (this  is  equivalent  to  X  =  0). 

r,  2 

Obviously,  this  approximation  is  too  inexact. 
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The  relaxation  mechanism  is  the  starting  point  of  the  data 


process.  Generally,  the  relaxation  mechanism  is  determined 
through  theoretical  analysis.  In  recent  years,  much  progress  [8] 
has  been  made  on  the  relaxation  mechanism  through  the  combination 
of  theories  and  experiments.  It  is  pointed  out  that  different 
mechanisms  could  possibly  exist  in  different  mixtures,  even  in  a 
the  same  system  under  a  different  temperature  condition. 

The  V-T  rate  which  is  derived  from  Schafer  formula  has  a 
certain  error.  It  is  not  suitable  to  use  the  Schafer  formula, 
especially  when  the  V-V  rate  is  not  much  greater  than  the  V-T 
rate.  The  formula  presented  in  this  paper,  however,  is 
appropriate  for  processing  this  type  of  data. 

Since  all  process  rates  are  calculated  from  a  certain  process 
data,  therefore,  if  the  data  used  as  the  base  are  not  accurate, 
the  calculated  rates  may  appear  as  large  errors  or  may  even  be 
nonsense  (such  as  a  negative  value).  Similar  situations  appear  in 
the  photosound  method.  Therefore,  one  has  to  be  careful  to  select 
the  rate  for  data  processing.  From  another  standpoint,  this 
provides  a  powerful  mean  to  examine  and  inspect  against  those 
existing  data  rates. 
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